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' §0- Introduction Here we consider difference equations of the form 

O 

. I \ P( z m Zn— l) /„ 1 N 

CN, z n+1 = g(z n ,z n -x) = r (0.1) 

^ : 

where g is a rational function of two variables, which we write as a quotient of polynomials 
; without common factor. Given starting values zq = x and z\ = y, equation (0.1) gives 
rise to an infinite sequence (z n ) n >o, as long as the denominator in (0.1) does not vanish. 
One of the basic questions here is to find the periodic sequences generated by (0.1), i.e., 
sequences for which Zj+ P = Zj for all j. We refer to the books [KL], [KM] and [GL] for 
expanded discussions of this question. Here we describe a method for giving an upper 
bound on the complex periodic sequences of period p. We note that equations (0.1) have 
been widely studied in the real domain, frequently for z n > 0; when we obtain an upper 
bound on the number of complex periodic points, this upper bound also applies to the 
periodic sequences which are real (or positive). 
Let us write (0.1) as a rational map 

O: F(x,y) = (y,g(x,y)) (0.2) 

| of the plane. It is most natural for us to consider F as a map of the complex plane C 2 . In 
fact, we will work with the extension of Fx '■ X — > X to a compactification X of C 2 . We 
let Ind(F) denote the set of points of indeterminacy of F. The set Ind(F) is finite and 
has the properties: (1) F is holomorphic on X — Ind(F), and (2) F cannot be extended 
continuously to any point of Ind(F). (In fact, F blows q up to a subvariety of X.) We 
will consider the set 

a - 

Fix' n = {z e C 2 : F J z eC 2 - Ind(F) V 1 < j < n, and f n z = z). 

Thus a point (x, y) belongs to Fix' n if and only if the corresponding sequence (z n ) is 
periodic for some period p which divides n. The mappings F that we consider here will 
have only isolated periodic points, and each such point q has a well-defined (algebraic) 
multiplicity n q , which at a regular point is the multiplicity of q as a solution of the pair of 
equations (x, y) — F n (x, y) = (0, 0). We set 

#Fix' n = 

qeFix' n 



Supported in part by the NSF. 



1 



denote the number of fixed points counted with multiplicity. The multiplicity is > 1, so 
an upper bound on j^Fix' n gives an upper bound on the number of points in Fix n . 
We will consider linear fractional recurrences of the form 

ao + a\z n -\ + a,2Z n 

Z n +1 = 7 - 7 —7 • 0.3 

b + biz n -i + b 2 z n 

For instance, the host-parasite model (see [KM, p. 181]), which is customarily presented 
as 

X n +1 = 



_ J^nyn 

Un+l ~. : n 

i + pyn 

may be re- written in the form (0.3) (see §2). Our first result is: 
Theorem 1. The periodic solutions of (0.3) satisfy 

#Fix' n <F n+ i + F n _i, (0.5) 

where F n denotes the nth Fibonacci number. For generic aj and bj, equality holds in (0.5). 
We may also consider the special case 

Z n +1 = r - • (0.6) 

b + z n -\ 

Let us define a Fibonacci-like sequence defined by (p n +3 = <p n + (p n +i and with starting 
values if i = 0, <f2 = 2, ^ps = 3. 

Theorem 2. The periodic points of (0.6) satisfy 

#Fix' n <<p n + 2. (0.7) 

Further, equality holds in (0.7) for generic a and b. 

The method of proof is an adaptation of the Lefschetz Fixed Point Theorem. In §1 
we discuss this method, which was used for a different family of rational maps in [BD]. In 
§2 we give the proofs of Theorems 1 and 2. And in §3 we show how it may be applied to 
two other examples. 

§1. Lefschetz Fixed Point Theorem Let us consider a polynomial map F : C 2 — > C 2 
which is given as 

F(x, y) = (f d (x, y),gd(x, y)) H , (1.1) 

where fd and are homogeneous polynomials of degree d, and the dots indicate terms of 
lower degree. We suppose that {fd = Qd = 0} = {(0,0)}. Such a polynomial map has a 
continuation to a holomorphic map F : P 2 — > P 2 , where P 2 denotes the projective plane, 
with points represented by homogeneous coordinates [xq : x\ : x 2 ]. Thus P 2 = C 2 U S , 



where we identify points of C 2 with points of P 2 by the map (x,y) h- > [1 : x : y], and 
Eq = {xq = 0} is the hyperplane at infinity. 

The cohomology of P 2 is #*(P 2 ; Z) = H°®H 2 ®H 4 = Z©Z©Z. For any rational map, 
the action of F* on H° is the identity and F* acts by multiplication by the topological 
degree <i t0 p on H A . Since F is holomorphic and is homogeneous of degree d, we have 
dtop = d 2 , and F* acts as multiplication by d on H 2 . Thus the total map on cohomology 
is given by the diagonal matrix 



F*\H*(P 2 ) = 




If d > 1, then the fixed points of F are isolated, and so the Lefschetz Fixed Point Theorem 
tells us that the number of fixed points is given by the trace of the mapping on cohomology: 

#Fix(F) = Tr(F*\H*(P 2 )) = l + d + d 2 . 

Here the expression #Fix(F) denotes the sum of the multiplicites of the fixed points 
in P 2 . If we wish to find Fixi(F), then we consider which of the fixed points are at 
infinity. We see that Eq is invariant, and the restriction F\T,q is given by the rational map 
[0 : x\ : x 2 ] !— > [0 : fd(xi,x 2 ) '■ gd{ x ii x 2)]- If we count the fixed points of F|Eq using the 
Lefschetz Fixed Point Theorem argument above, we find #Fix(F\Y, ) = 1 + d. Thus we 
can count the fixed points in C 2 as #Fix[(F) = d 2 . 

Now we observe that F n is a mapping satisfying (1.1), and the degree of F n is d n . 
Applying our formula, we have #Fix' n = d 2n . 

Let us next consider a mapping that arises in the so-called SI epidemic model (see 
[KM, page 186]): 

g:(S,I)^(S + aSI,I-aSI). (1.2) 
We note that a related model, of the form 

(S, I)^(S + aSI + XI, (1 - X)I - aSI) (1.3) 

can be treated in the same way. We may extend g to projective space by setting S = x\/xq 
and I = X2/X0. This gives 

g[x : x\ : x 2 ] = [xq : Xi(x - ax 2 ) : x 2 (xo - axi)]. 

We see that g defines a holomorphic map of P 2 except at the indeterminacy locus 

/ne*(<7) = {ei = [0:l:0],e 2 = [0:0:1]}. 

We see that both points e\ and e 2 blow up to the hyperplane at infinity E . The hyperplane 
at infinity is exceptional in the sense that 

E -{ei,e 2 } ^[0:1: -1]. 
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Further, we note that [0:1: —1] G E and is a fixed point of g. 

We observe that the algebraic degree of g is 2, and the topological degree is also 2. 
Thus we have 



g*\H*(P 2 ) = | g*\H 2 

degtop 

Now the trace of g* is equal to the intersection number between the diagonal A C P 2 x P 2 
and (7 _1 A (see [Fu]). In the case of our mapping, the intersection points are isolated, and 
we have 

Tr(g*) = 1 + 2 + 2= £ N . 

qeAng-^-A 

The points of A n g _1 A — Ind(g) are exactly the fixed points of g. For a point 
g G (A fig -1 A) fl Ind(g) there are two possibilities: (1) (/ blows up to a curve 7 E5 q; or (2) 
there is an exceptional curve E which contains q and g(E — Ind(g)) = q; (or both). We 
have seen above that e\ and e2 are in case (1): they are indeterminate points which blow 
up to Eq, and e±, G Eq. We note that while [0:1: —1] ^ Ind(g) it is in case (2), since 
it is the image of the exceptional curve Eq. Further, since Eq is mapped to itself, we see 
that there can be no orbit which intersects both C 2 and Eq. Thus there are no periodic 
points whose orbits pass through infinity. This gives: 

#Fix' 1 (g) = Tr(g*\H*) -3 = 2. 

We observe that applying this argument to (1.2) and (1.3) is the same, since we are only 
concerned with the behavior of points at infinity, and this is not influenced by the lower 
degree terms of g. 

In order to discuss Fix' n , we consider the iterate g n . Eo is the only exceptional curve 
for g, and Eo is taken to a fixed point. Thus there is no exceptional curve that is mapped 
to Ind(g) by an iterate of g. This is the criterion given in [FS] for a map of P 2 to satisfy 
deg(g n ) = (deg(g)) n = 2 n . Thus we have 



\H* = 



1 

2 n 

2 n 



and so the intersection number of A fl g~ n A is 1 + 2 • 2 n . Now the 3 points e±, ei and 
[0 : 1 : — 1] exhaust all the indeterminacy points and exceptional images, so there are no 
new "fake periodic points" for higher period. This gives 

#F<(<7) = 2.2™-2. 



§2. Linear Fractional Recurrences We wish now to show how to apply the ideas of 
§1 to the maps (0.3) and (0.6). Let us start by rewriting the recurrence (0.3) as a planar 
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map of the form (0.2). Rewriting this mapping in terms of homogeneous coordinates 
[xq : X\ : X2] = [1 : x : y], we have 



where a = (0,0,0,1,0,2) and (3 = (&0;&i;&2)- F is invertible, so the topological degree is 1, 
and we have 



The points of indeterminacy are Ind(F) = {ei,po,p-y}, where po = [0 : — (3% '■ Pi] and 
P-y = [P1&2 — '■ —f3o&2 + «o/32 : Oiifio — olq(3\\. See [BK, §1] for details. For generic 
values of a and (3, e\ G A n / _1 A has multiplicity 2. First, ei is a fixed point because 
ei G So, which is an exceptional curve, and £0 f— ¥ e i 5 an d second it satisfies the property 
(2) above, because e± blows up to the curve and ei G Sb- Taking the intersection 
number and subtracting the multiplicity of e\, we have 



We note that for non-generic values of a and (3 this becomes #Fix\{F) < 2. 

We now want to treat higher iterates of F, but deg(F n ) 7^ (deg(F)) n . This is seen 
because, Eo ^ ei G Ind(F), i.e. the forward orbit of an exceptional curve enters the 
indeterminacy locus. We deal with this situation by using the approach of Diller and 
Favre [DF], in which we replace P 2 by a space on which the passage to cohomology is 
consistent with iteration. In this case, we work with the space Y which is obtained by 
blowing up the point e±; further details are given in [BK, §2,3]. H 2 (Y;Z) is generated 
by the class of a hyperplane H and the blowup fiber E\. With respect to this basis, the 
induced rational map F Y : Y — > Y acts on H 2 (Y; Z) according to 



The map Fy has two exceptional curves, which are denoted and E 7 , and Ind(Fy) 
consists of two points. For generic a and (3, the orbits of and E 7 do not meet the 
indeterminacy locus, so the map Fy is regular on H 2 , which means that (Fy)* = (F Y ) n . 
Using this, as well as the defining properties of the Fibonacci numbers F n , we have 



F[xq : X\ : X2] = [xq(3 ■ x : X2(3 ■ x : XqCx ■ x] 




#Fzx;(F) = Tr(F*\H*) -2 = 2. 





Thus we obtain 



#Fix' n = F n+2 - F 7 



n 



2 — F n+ i + F n _i. 



This proves Theorem 1. 
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Another formulation. Let us show how to rewrite the mapping (0.4) as a (linear 
fractional) difference equation (0.3) with ao = a\ = 0, a 2 = a, bo = 1, b\ = 7, and 
b 2 = —13 /a. Applying (0.4) we have 

a ix n + ^l3x n y n ^I3x n y n jx n 7 

IXn - PVn+1 = r— ~ = - = 

1 + py n 1 + /?y n 1 + PVn ot 

Thus 

7 1 

J/n+1 = -j(^n X n+1 ), (2.1) 

which means that the sequence {(x n ,y n )} in (0.4) is determined by the sequence {x n }. 
From the first equation in (0.4) and (2.1) we have 

Xn+l 



1 + 7x n _i - ^xj 
so the sequence {x n } satisfies (0.3) with the given values of cij and bj. 

Proof of Theorem 2. Now let us consider the map (0.6), which is the mapping (0.3), 
but restricted to the (non-generic) case a = (a, 0,1) and (3 = (6,1,0). There are three 
points of indeterminacy Ind(F) = {ei,e 2 ,p}, and there are three exceptional curves Eo, 
Tip and E a . We have Tp 1— > e 2 G Ind(F) and Eo 1— > ei G Ind(F). We let X denote the 
space obtained from P 2 by blowing up ei and e 2 . We let Fx ■ X — > X denote the induced 
rational map. Fx has one indeterminate point Ind{F x ) = {p}. The line E a C I is 
exceptional, and E a 1— > q 1— > -FxQ 1— > Fj^-g 1— > • • •. For generic values of a and 6, the orbit of 
the exceptional line does not land on Ind(Fx) = {p}- In this case we have (F x ) n = (F x )*. 

If H denotes the class of a general hyperplane in X , and if E\ and E2 are the excep- 
tional blowup fibers, then {H, Ei,E 2 } gives an ordered basis for H 2 (X; Z). In [BK, §4] we 
saw that 

/ 2 1 1 
F x = -1 -1 

V-i -1 -1 

with respect to this basis. 

The characteristic polynomial of A := F x is x s — x — 1 so we have A 3 — A — I = 0. It 
follows that each entry a™ ■ of A n satisfies 

Now let <^ n be the sequence satisfying the linear recurrence tp n +3 = <fn+i + and which 
has starting values cp = 3 (= Tr(J)), ^ = (= Tr(A)), and ^ 2 = 2 (= Tr(A 2 )). It 
follows that y> n = Tr(A n ) for all n, which proves Theorem 2. 

§3. Two Examples In this section we will consider two maps that arise as biological 
models. The first (see [KM, p. 174]) is 

g: (x n+ i,y n+ i)= ax n H ■ ■ ,py n + 



a + aix n + a 2 y n b + b\x n + b 2 y n , 
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For this map, we determine the topological degree as well as the degree as a homogeneous 
polynomial. Then we show that deg(g n ) = (deg(g)) n and count the spurious fixed points. 

Finding the number of preimages of a generic point, we see that d top = 4. Checking for 
indeterminate points in C 2 , we see that we have fractions of the form jj at p a := (0, —0,0/0,2) 
and pb := (— bo/bi, 0). The point p a blows up to a horizontal line, which contains e\ £ P 2 ; 
and pb blows up to a vertical line, which contains e 2 £ P 2 . 

If we write g in homogeneous coordinates, we find the mapping 

g[x : x\ : x 2 ] = [x £ a (x)£ b (x) : x 1 (. b {x)m 1 {x) : x 2 £ a (x) m 2 {x)\ 

where £ a {x) = oqXq + a\X\ + a 2 x 2 , lb is similar, and m\ and m 2 are linear functions. Thus 
g has algebraic degree 3. Let us set L a = {l a = 0} and Lb = {lb = 0}. The points of 
indeterminacy are then Ind(g) = {p a ,Pb,Qa,Q.b, r }, where q a = L a fl £0 = [0 : — a 2 : ai], 
q b = Lb fl E = [0 : —b 2 : 61], and r = L a fl We observe, also, that p a £ L a and 
p b £ Lb. The Jacobian determinant of g is given by i a ibR{x), where R(x) is quartic. In 
fact, {R = 0} is a curve of branch points. Thus the set of exceptional curves consists of 
L a l— ei and L5 1— > e2- The points ei and e2 are regular (not indeterminate) and are fixed. 
Since the orbits of exceptional curves never enter the indeterminacy locus, it follows that 

/I 

deg(g n ) = (deg(g)) n = 3 n . Thus we have g n *\H*(V 2 ) = 3 n 

Now we count the points that are not in Fix' n , but which will contribute to the 
intersection number of A fl / _1 A. Namely, we have the (regular) fixed points e\ and e 2 
at infinity. Further, we have p a , which blows up to L a 3 p a , and p b , which blows up to 
Lb 3 Pb, for a total count of 4. For generic parameters, we have q a 7^ qb, etc, so taking the 
trace and subtracting the excess multiplicity, we have 

#Fix' n = 3 n + 4 n - 3. 

We apply a similar analysis to 

/ axy + x + y cxy + x + y 



f{x,y) = 



\ b+x+y ' d+x+y 



(see [KM, p. 172]). This map differs from the previous one because it needs to be regular- 
ized. 

First, we note that the topological degree of / is 2. If we rewrite this mapping in 
terms of homogeneous coordinates [xq : x\ : x 2 ] = [1 : x : y] on projective space P 2 then 
we have 

f[x \x\ : x 2 ] = [x (bx + x\ + x 2 )(dx + x\ + x 2 ) : 

(dxo + x\ + x 2 ){ax\x 2 + x$x\ + xqx 2 ) : (bx + x\ + x 2 )(cxix 2 + x xi + x x 2 )]. 

To find the exceptional curves, we look at the Jacobian, which is 

x (xi - x 2 )(bx + xi + x 2 ){dxQ + x\ + x 2 )Q(x) 
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for some quadratic polynomial Q(x). The line {x\ = x 2 } consists of branch points; all of 
the other critical curves are exceptional: 

{xo = 0} I— > [0 : a : c] E E , {bxo + x\ + X2 = 0} i— > ei, 

{<ixo + cci + ^2 = 0} I— > e 2 , {<5(a^) = 0} i— > [ac(6 — <i) : a(a — c)d : 6(a — c)c]. 

We observe that [0 : a : c] G Eq is a fixed point of / (at infinity). The three points of 
indeterminacy are obtained by finding the points of intersection of exceptional curves: 

/nd(/) = {[0:l:-l],ei,e 2 }. 

We see that [0:1: —1] blows up to E 3 [0 : 1 : —1], so that this point will also be part 
of our intersection count when we are using the Lefschetz index. 

Now let 7T : Z — > P 2 denote the complex manifold obtained by blowing up e\ and e 2 . 
If H is the class of a generic hyperplane, and if E\ and E 2 are the exceptional (blow-up) 
fibers, then {H, E±, E 2 } is a basis for H 2 (Z). Since e± and e 2 both blow up to lines, we 
find that f%H = 3H — E\ —E 2 . Further, the preimages of e\ and e 2 are lines which contain 
neither e\ nor e 2 . Thus f%E 1 = f* z E 2 = H. 

It follows that 



r z \H*{z) 



f 1 



\ 



i 






1 






\ 



2/ 



Now we consider the behavior of fz '■ Z — > Z. The only exceptional hypersurfaces 
are {xq = 0}, which maps to a fixed point, and {Q(x) = 0}, and there is only one point 
of indeterminacy: Ind(fz) = [0:1: — 1]. For generic parameters, the orbit of {Q = 0} 
will be disjoint from Ind(fz)- Thus, for generic parameters, we have (f z ) n = (fz)* f° r 
n = 1, 2, 3, . . . The number of spurious periodic points for fz is still 2, and so we subtract 
them from the index to obtain 



#Fix' n = Tr(f* n ) - 2. 
The characteristic polynomial of f z \H 2 (Z) is x(x 2 — 3x + 2), so r n 



Tr(f z n \H 2 



satisfies r n + 2 = 3r n +i — 2r n , with initial condtions t\ = 
and for generic parameters we have #Fix' n = r n + 2 n — 1 
this gives an upper bound for #Fix' n . 



) 

3 and r 2 = 5. Thus r n = 2 n + 1, 
= 2 n+1 . For arbitrary parameters, 
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